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Superconducting High Pressure Phases Composed of Hydrogen and Iodine 
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Evolutionary structure searches predict three new phases of iodine polyhydrides stable under pressure. Insu¬ 
lating P1-H5I, consisting of zigzag chains of (HI)^+ and molecules, is stable between 30-90 GPa. Cmcm- 
H 2 I and P 6 /mmm-H 4 l are found on the 100, 150 and 200 GPa convex hulls. These two phases are good 
metals, even at 1 atm, because they consist of monoatomic lattices of iodine. At 100 GPa the Tc of H2I and H4I 
are estimated to be 7.8 and 17.5 K, respectively. The increase in Tc relative to elemental iodine results from a 
larger cciog from the light mass of hydrogen, and an enhanced A from modes containing H/I and H/H vibrations. 

PACS numbers: 74.70.-b, 74.62/Fj, 62.50.-p, 63.20.dk, 61.50.Ks 


We have only begun to uncover the affects of pressure on 
the stoichiometries, electronic structures and emergent behav¬ 
ior of condensed matter fDID- A number of phases with 
combinations that appear strange from a 1 atm perspective 
have been theoretically predicted and experimentally synthe¬ 
sized. These include hydrides with unique combinations such 
as LiH 2 and LiHe O [H, as well as H^S phases 111 , some 
of which may have a superconducting transition temperature, 
Tc = 190 K at 150 GPa |6H9l, that surpasses those of the 
oxide-based high-temperature superconductors. Other pre¬ 
dicted phases include sub- or polyhydrides containing an al¬ 
kali or alkaline earth metal EiiMia and hydrides contain¬ 
ing a Group 14 element |[T9tiT7l . Theory has suggested that 
some of these hydrides may be superconducting at high tem¬ 
peratures CD, or contain unusual structural motifs such as the 
linear symmetric H 3 anion - the simplest example of a three- 
center-four-electron bond m. 

Another example of unique stoichiometries that could po¬ 
tentially be synthesized under pressure are H 2 CI and H 5 CI 
l28l . They were predicted to become stable with respect 
to decomposition into H 2 and HCl at 60 GPa, and remain 
so until at least 300 GPa. All of the H^Cl phases con¬ 
tained zigzag [HCl] 00 chains reminiscent of those present 
within HCl. The three-center-two-electron (3c-2e) H^ mo¬ 
tif, which approached the ideal equilateral triangle configura¬ 
tion by 300 GPa, was incorporated within H 5 CI. How would 
decreasing the electronegativity of the dopant atom affect the 
stable crystal lattices and their electronic structures? To an¬ 
swer this question we theoretically studied the iodine polyhy¬ 
drides, H^I. A further motivation for our work: it has been 
proposed that pressure can be used to synthesize alloys of 
metallic hydrogen with the general formula (H 2 )i-a;HIa^ that 
are potentially superconducting |29l . Moreover, HI undergoes 
a pressure induced insulator to metal transition below 50 GPa 
I 29 I . and iodine becomes metallic at 16 GPa 1^ and super¬ 
conducting below 1.2 K at 28 GPa 1311 . 

The isomorphic low temperature phases of HF, HCl 
and HBr contain planar zigzag chains of hydrogen-bonded 
molecules held together by weak van der Waals (vdW) forces. 
HI assumes a planar distorted hydrogen-bonded diamond lat¬ 
tice, but its detailed structure is not known l32l . Density 
functional theory (DFT) calculations have been undertaken to 


study the structural evolution of the heavier hydrogen halides 
up to 200 GPa |33]| . But, because the phases of HI under pres¬ 
sure are unknown we proceeded to find the global minimum 
structures from 0-200 GPa in 50 GPa intervals using the evo¬ 
lutionary algorithm (FA) XtalOpt |34l coupled with PBE- 
DFT l35]| calculations carried out with VASP |36l. More in¬ 
formation about the computational details is provided in the 
Supplemental Material, SI ITtI . The most stable structures 
that emerged in our FA searches were comprised of 2D segre¬ 
gated layers of H 2 and iodine, implying that HI is metastable 
with respect to these constituents. In addition, we optimized 
geometries where HI assumed the experimentally determined 
and theoretically predicted low-temperature structures of the 
lighter hydrogen halides |33l. All of these were less stable 
than the H 2 /I segregated phases. These results suggest that 
the impurities and side-products observed upon compression 
of HI above 70 GPa |38]| may have been due to decompo¬ 
sition of the halide. A previous theoretical study found that 
HBr was unstable with respect to decomposition into H 2 and 
Br 2 between 120-150 GPa |3^ . At these pressures bromine 
is monoatomic, as is iodine, and it should not be unexpected 
that HI behaves similarly but at lower pressures. 

FA runs were carried out to find the most stable crystal lat¬ 
tices of H^I (n = 2-9) at 50, 100, 150 and 200 GPa. The struc¬ 
tural motifs present in the hydrogen-rich iodine phases that 
emerged from our searches fall into one of two categories: 
1 ( 5 +, molecules as found in H 5 I, or 1^^+ and 

H 2 “ molecules as found in H 2 I and H 4 I. In Ref. 1^ it was 
shown that density functionals approximating the effects of 
vdW interactions, and zero-point-energy contributions to the 
enthalpies both altered the transition pressures between the 
most stable H^Cl phases slightly, but they did not affect the 
identify of the phases nor their stability. For this reason we 
also expect them to have a negligible effect on our results. 

Fig. shows the AHf of the predicted phases with respect 
to the elemental solids. The convex hull is defined as the set of 
line segments below which no other points lie, and the phases 
whose AHf comprise the hull are thermodynamically stable 
with respect to decomposition into other polyhydrides and/or 
solid H 2 and iodine. At 50 GPa H 5 I is the only species on the 
hull. This stoichiometry becomes stable with respect to de¬ 
composition into H 2 and I at ^30 GPa and it has the most neg- 


2 



FIG. 1: (a) The enthalpy of formation, AiT/, for the reaction 
(|H2)n+I ^ Hnl VS. the H 2 composition as a function of pressure. 
For H 2 we used the enthalpy of the structures from Ref. |[^ and for 
iodine the structure sequence Cmca Immm lA/mmm 
Fm3m ll40lB3l . The computed transition pressures of 18, 42, and 
55 GPa are in good agreement with the experimental ones, 21, 43, 
55 GPa (neglecting the incommensurately modulated phase 1441 ). 
Fm3m iodine persists until at least 275 GPa 103). (b) AHf vs. 
pressure of select H^I phases. becomes negative for 

P > 30 GPa. Above 90 GPa H 2 I has the most negative AHf. 


ative AHf of any of the structures examined until 90 GPa. At 
this pressure the H 2 I stoichiometry becomes the lowest point 
on the hull and remains so until at least 200 GPa. H 4 I also 
comprises the 100, 150 and 200 GPa hulls. At 200 GPa the 
AHf of all of the phases continues to decrease. Even though 
some of the stoichiometries predicted to be stable in our cal¬ 
culations were also found in the hydrogen rich H/Cl phase 
diagram under pressure ||28|, the structures adopted and their 
properties are distinct, highlighting the importance of the elec¬ 
tronegativity (2.66 for I and 3.16 for Cl on the Pauling scale) 
and radius (1.98 A for I and 1 .75 A for Cl) of the halogen atom 
on the emerging structures and their properties. Importantly, 
whereas the H^Cl structures were not superconducting 1281 . 
we predict that electricity may pass without resistance through 
H 2 I and H 4 I below ~7.8 and 17.5 K at 100 GPa, respectively. 

The most stable H 5 I phase possessed Pi symmetry up to 
90 GPa, and it assumed Pmma symmetry at higher pressures. 
Above 150 GPa another phase, containing only molecular hy¬ 
drogen, became preferred, but it did not fall on the convex 
hull. Phonon calculations at 50 GPa verified the dynamic sta¬ 
bility of P 1 -H 5 I. Just like the H^Cl phases predicted in Ref. 
flSl , H 5 I contained zigzag [HX]oo (X=C1, I) chains that re¬ 
sembled those found in HCl, HBr, and HI at low temperatures 
Ea. These chains lay parallel to the crystallographic 6 -axis 
in P 1 -H 5 I, as illustrated in Fig.[^ At 50 GPa the H-I distances 
along the chains were nearly equivalent, 1.795 and 1.788 A, 
indicative of pressure induced multi-center bonding, which 
has also been observed in HX (X=F, Cl, Br) 13^1 and H 2 CI 
(281. The H-I-H and I-H-I angles measured 85.7-85.8° and 



FIG. 2: (a) Supercell of P1-H5I at 50 GPa. Hydrogen atoms are 
white, iodines are purple. The zigzag HI chains are a common motif 
in compressed hydrogen halides, (b) Contour plot of the ELF map of 
Pi -H5I shown in a plane parallel to the 6-axis, (c) Total and projected 
DOS of P1-H5I at 50 GPa. Ep has been set to zero. 


179.4-179.8°. The chains arranged to form sheets that were 
surrounded by layers of H 2 units in which the H-H bond mea¬ 
sured 0.747 A, which is nearly identical to what we calculate 
for the isolated gas phase molecule at 1 atm, 0.750 A. Actu¬ 
ally, six H 2 molecules surround half of the hydrogen atoms 
in the [HI] 00 chains in a “side-on” fashion. In our previous 
studies we found a number of structures where H 2 encircled 
the electropositive metal in this same manner, for example in 
MgHi 2 and MgHie ca. A Bader charge analysis revealed 
that the charge on the hydrogen and iodine atoms in the H-I 
chains are +0.03 and +0.12, respectively. This can be com¬ 
pared to - 0.01 and + 0.01 for the hydrogen and iodine atoms in 
a hypothetical HI phase at 50 GPa whose Cmc2i symmetry 
was chosen because HE, HCl and HBr assume this structure 
at ambient pressure. All of the hydrogen atoms in the H 2 
molecules were found to bear a slight negative charge such 
that the average is (H“^-^^) 2 . The distance between the hy¬ 
drogen atoms comprising H 2 and the hydrogen atoms in the 
H-I chains ranged from 1.950-2.022 A and the closest iodine 
atom was at least 2.273 A away. This suggests that the H 2 
molecules do not interact strongly with any of the atoms com¬ 
prising the [HI ]00 zigzag chains. Indeed, the H 2 vibron fre¬ 
quency ranged from ^4150-4250 cm“^, which compares well 
with the 4161 cm“^ observed for gaseous H 2 at 1 atm na. 
To probe this further, we calculated the electron localization 
function (ELE) for P1-H5I at 50 GPa, see Eig.j^b). The ELE 
reveals a high tendency of electron pairing within the H^” 
units and along the nearly symmetric [HI ]00 chains. In fact, 
the ELE along the chains is reminiscent of the one calculated 
for symmetric HBr at 30 GPa (3^ . 

CC-H 5 CI was found to be a stable structure in the hydro¬ 
gen rich H/Cl phase diagram between 100-300 GPa (28l. By 
300 GPa H^, the simplest example of a 3c-2e bond, emerged 
in this phase. The bond lengths (^0.89 A) and angles (60°) in 
the solid approached those in the gas phase species, and sig¬ 
natures of H^ vibrations were observed in the phonon DOS 
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between 1800-3400 cm“^. We do not see any evidence for 
the formation of within H 5 I up to at least 200 GPa. 

At 50 GPa P 1 -H 5 I is insulating with a PBE band gap of 
-0.5 eV, see Fig. [gb). The I s-states overlap primarily with 
the 5 -states, whereas I p-character is found in a region 
containing both and contributions. Metalization 
occurs as a result of pressure induced broadening of the va¬ 
lence band, primarily of iodine p-character, and the conduc¬ 
tion band, which exhibits mainly 5 -character, by 65 GPa 
within PBE. But, because GGA functionals underestimate the 
pressure associated with band gap closure it is likely that this 
phase would metalize at a pressure outside the stability range 
of the H 5 I stoichiometry. For comparison, HCl and HBr were 
computed to be metallic by 130 and 80 GPa, and the bands at 
Ep were primarily of halogen p-character l(33]| . 

Between 90-200 GPa a H 2 I phase with Cmcm symmetry 
had the most negative AP/, and phonon calculations revealed 
that it was dynamically stable at 100 GPa. This structure was 
found to be the most stable H 2 I arrangement in evolutionary 
searches performed at 50-200 GPa. Two rows of H 2 molecules 
arranged in a zigzag fashion run through channels formed by 
the iodine host within this structure, see Fig.j^a). At 100 GPa 
the hydrogen molecules measured 0.796 A and the distance 
between them 1.399 A, as compared with 0.732 and 1.495 A 
in solid H 2 at this pressure. The I-I contacts measured be¬ 
tween 2.896-2.904 A, which is very close to the 2.899-2.920 A 
we compute within Fm3m iodine at 100 GPa. The shortest 
H-I contacts were 2.097 A, suggesting little interaction be¬ 
tween the constituents of the two lattices. Cmcm- 112 ^ bears 
no resemblance to the H 2 CI phase that was predicted to be 
stable by 100 GPa |[28l. High ELF values were calculated 
between the hydrogens comprising the molecules. Simi¬ 
lar to compressed solid iodine, a large ELF between the iodine 
atoms was not observed. A Bader analysis showed that charge 
is transferred from iodine to hydrogen, (H “^-^^)2 at 100 
GPa. By 200 GPa the H 2 bonds within Cmcm-1121 stretch 
slightly to 0.817 A and the distance between them shrinks to 
1.198 A. This suggests that if Cmcm-ll 2 l is stable at higher 
pressures, these chains may polymerize resulting in a 1-D 
hydrogenic motif, as has been predicted for RSm-SrHQ by 
250 GPa (nl. 

A P6lmmm symmetry H 4 I structure also comprised the 
100 GPa convex hull, and it was found to be the most stable 
H 4 I phase between 90-200 GPa. Phonon calculations revealed 
dynamic stability at 100 and 150 GPa. This phase is com¬ 
posed of 1-D chains of iodine atoms with an I-I distance of 
2.812 A and H 2 molecules measuring 0.799 A, both lying par¬ 
allel to the c-axis as illustrated in Fig.j^b). The intermolecular 
distances between the dihydrogen molecules (1.875-2.013 A) 
and between H-I (—2.13 A) are too long for bonding interac¬ 
tions. Similar to H 2 I, the iodine donates electrons (the Bader 
charges are (H“^-^^) 4 l^-^^), and the ELF indicates high elec¬ 
tron localization only within the dihydrogen molecules. At 
200 GPa the H 2 units are not yet close enough to polymerize. 

The total and projected DOS plots for Cmcm-ll 2 l and 
P 6 /mmm-H 4 l at 100 GPa resemble each other, see Fig. 
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FIG. 3: Supercells and the total and projected DOS plots of (a) 
Cmcm-ll 2 l and (b) PG/mmm-H^l at 100 GPa. 


which is not surprising considering their structural similari¬ 
ties. Both phases are good metals with the DOS at Ep being 
0.061 and 0.049 eV“^ states/(valence electron), respectively. 
The DOS can be split up into 3 regions: <-13 eV, -13 to -5 eV 
and above -5 eV it is primarily of I 5, H 2 5 and I p-character, 
respectively, but the H 2 5- states have a large degree of overlap 
with both iodine 5 and p. The metallicity is primarily due to 
the iodine p-states, but for P >100 GPa the hydrogenic states 
also provide a substantial contribution to the DOS at Ep. 

Pressure causes the intramolecular I-I bond in Cmca-l 2 
to stretch, and the intermolecular distances to decrease. At 
— 16 GPa this phase undergoes an insulator to metal transi¬ 
tion 1 ^ . which is superseded by an incommensurate phase 
at 23 GPa and the formation of a monoatomic phase within 
which the I-I distances are “equalized” by 30 GPa il. Ex¬ 
periments yielded a bond length of 2.75 A at ambient condi- 
tions, whereas the nearest neighbor bonds in the fully disso¬ 
ciated monoatomic crystal measured 2.89 A 1441 . In our cal¬ 
culations bond equalization occurs near 35 GPa with I-I dis¬ 
tances of 2.98 A. Interestingly, we find that Cmcm-H 2 l and 
P 6 /mmm-H 4 l are already metallic at 1 atm, and the DOS 
ai Ep exhibits primarily iodine p-character. The metallicity 
is a result of the monoatomic iodine present in these phases 
at ambient conditions, wherein the nearest-neighbor distances 
are found along 1 -D chains measuring 2.93 A. 

Experimental measurements yield a Tc of 1.2 K at 28 GPa 
for iodine ED Theory has shown that under hydrostatic con¬ 
ditions Tc decreases with increasing pressure, falling to < 1 K 
above 100 GPa (461 1121 • Our calculations for the Em3m 
phase of iodine using the Alien-Dynes modified McMillan 
equation (481 yield Tc values that are in agreement with these 
studies, see TableThey further confirm that the logarithmic 
average frequency, cciog, which is low due to the heavy mass of 
iodine, increases slightly with pressure ED- But the decrease 
in the electron-phonon coupling. A, with increasing pressure 
dominates, leading to a drop in Tc at higher pressures. 

What effect would incorporation of H 2 within lattices con- 
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TABLE I: Electron-phonon coupling parameter (A), logarithmic av¬ 
erage of phonon frequencies (ct;iog) and estimated superconducting 
critical temperature (Tc) for values of the Coulomb pseudopotential 
(//*) of 0.1 and 0.13 for Cmc?Ti-H 2 l, PG/mmm-K^l and Fm3m 
iodine at various pressures. 


System Pressure (GPa) 

A 

c^iog (K) Jf 

= ° i (K) Jf 

= 0.13 

H 2 I 

100 

0.51 

608.7 

7.8 

4.8 

H 4 I 

100 

0.58 

854.3 

17.5 

12.0 

H 4 I 

150 

0.56 

1097.4 

20.4 

13.7 

I 

100 

0.31 

251.6 

0.16 

0.03 

I 

150 

0.25 

304.5 

0.02 

0.00 


sisting of monoatomic iodine, in the geometries adopted by 
H 2 I and H 4 I, have on cjiog, A and Tc? At 100 GPa we find 
that cjiog becomes larger as the H:I ratio increases, see Table 
|l| because of the light mass of hydrogen. In addition, the A 
of H 4 I is slightly larger than that of H 2 I, and these are both 
significantly larger than that of pure iodine. The most signifi¬ 
cant contribution towards A in both of the iodine polyhydrides 
stems from vibrational modes between 450-1300 cm“^ that 
involve both hydrogen and iodine motions, with the hydrogen 
vibron also contributing substantially towards A, see Fig. 
Vibrations that involve only the iodine atoms do not contribute 
much to the total electron phonon coupling. The computed Tc 
of the iodine polyhydrides for /i* = 0.1, 7.8 and 17.5 K for 
H 2 I and H 4 I at 100 GPa, is significantly larger than that of 
monoatomic iodine because H 2 increases both cciog and A. 

We also examined the effect of pressure on the Tc of H 4 I. At 
150 GPa the Tc was computed to increase to 20.4 K as a result 
of the larger cciog, which is due to the larger frequency range 
of the iodine and mixed iodine-hydrogen vibrational modes. 
In contrast, the frequency associated with the H 2 vibron de¬ 
creased slightly with increasing pressure as a result of bond 
weakening within the hydrogen molecule, which measured 
0.794 and 0.806 A at 100 and 150 GPa, respectively. A Bader 
analysis revealed that the softening of the hydrogen vibron can 
be attributed at least in part to the increased electron donation 
from iodine to at higher pressures, which lengthens this 
bond. The slight decrease in the electron-phonon coupling 
with increasing pressure stems primarily from the smaller val¬ 
ues of A associated with modes arising from motions involv¬ 
ing iodine atoms. In contrast, the A associated with the H 2 
stretch increased slightly under pressure. 

In conclusion, evolutionary searches have identified three 
hitherto unknown iodine polyhydride phases with unique stoi¬ 
chiometries that become thermodynamically and dynamically 
stable above 30 GPa. P 1 -H 5 I, composed of zigzag [(HI)‘^+]oo 
chains and molecules, is found on the 50 GPa convex 
hull and has the lowest AHf of any phase examined until 
90 GPa. Metalization occurs as a result of pressure induced 
broadening of the iodine p and hydrogen s-states, at pressures 
likely outside of the stability range of this structure. The DOS 
at Ef suggests that the Tc of this phase is likely to be low. 
Cmcm-H 2 l and P 6 /mmm-H 4 l both lie on the 100, 150 and 
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PIG. 4: Phonon band structure, the Eliashberg spectral function, 
a‘^F{uj), and the electron-phonon integral, A(cc;), for FGlmmm- 
H4I and Cmcm-ll 2 l at 100 GPa. The insets provide the contribution 
to A (absolute value and percent) from vibrational modes that are 
comprised of iodine vibrations (I-I), motions of hydrogen and iodine 
atoms (H-I), and the H2 vibron (H-H). 


200 GPa convex hulls. They are comprised of monoatomic 
iodine lattices that render these phases metallic by 1 atm, and 
H 2 molecules. A Bader analysis suggests a slight transfer of 
charge from iodine to H 2 . The metallicity persists as H 2 I 
and H 4 I become dynamically and thermodynamically stable 
phases. The DOS at Pi? is primarily due to the iodine p-states, 
with an admixture of hydrogen s. At 100 GPa the Tc of H 2 I 
and H 4 I is estimated as being 7.8 and 17.5 K, respectively. 
The presence of H 2 significantly enhances Tc in comparison 
to pure iodine because the light atomic mass of hydrogen in¬ 
creases cciog, and because vibrations involving H and I atoms, 
as well as the H 2 vibron, contribute substantially to A. As 
the quest for high-temperature superconductivity in hydrides 
gains momentum, we look forward to the eventual synthesis 
of these, and other intriguing phases. 
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